Understanding the Effect of Atmospheric Density on the Cosmic Ray Flux Variations 

at the Earth Surface 
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We report in this letter for the first time the numerical simulations of muon and neutron flux 
variations at the surface of the earth with varying air densities in the troposphere and stratosphere. 
The simulated neutron and muon flux variations are in very good agreement with the measured 
neutron flux variation in Oulu and the muon flux variation in Atlanta. We conclude from this study 
that the stratosphere air density variation dominates the effects on the muon flux changes while the 
density variation in troposphere mainly influences the neutron count variation. These results pave a 
new path for systematically studying the global temperature evolution using worldwide cosmic ray 
data. 

PACS numbers: 92.60.hv, 92.70.Kb, 92.60.Wc, *91.62.Xy 
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Sporadic weather patterns occur so frequent in recent 
years. According to NASA, the year 2012 was the ninth 
warmest year since 1880. The temperature at the sur- 
face of the earth has increased by ~0.6°C during the last 
century [![. Currently, there are two competing inter- 
pretations of the causes for the global warming: natural 
forcing such as solar influences on the climate, and influ- 
ences of anthropogenic activities. The studies done by 
Svensmark Q and others [1, H[, on the other hand, in- 
dicate that the galactic cosmic rays (GCR) may play a 
significant role on the temperature variation of the earth. 

Over the past decades, quite a few studies reported on 
the correlations between the earth weather and cosmic 
ray flux 0, It is showed that an 11-year average 

of northern hemispheric land and marine temperature is 
more closely correlated with the measured variation in 
cosmic ray flux [2]. Furthermore, there exists a causal 
relationship between GCR and low cloud coverage (< 3.2 
km) which suggests that the formation of the low clouds 
is influenced by the atmospheric ionization produced by 
GCR @. Clouds play a significant role on the earth's 
radiation budget by reflecting the incoming short waves 
and trapping the outgoing long waves. This causal rela- 
tionship reflects that the variation in GCR may indirectly 
influence the earth's temperature. 

While the true impact of cosmic rays on the earth cli- 
mate change is currently far from conclusive, continued 
efforts of long-term monitoring of cosmic ray flux varia- 
tions are imperative. This study also requires a quantita- 
tive understanding of the influence of atmosphere air den- 
sity fluctuations to the cosmic ray flux. Figure Q] shows 
the neutron daily counts percentage variation (blue) and 
pressure (red) measured at the Sodankyla Geophysical 
Observatory in Oulu, Finland in 2012. The Sodankyla 
Geophysical Observatory has been recording the neutron 
counts since 1964. At Georgia State University, a sea- 
sonal muon flux variation is observed as shown in Fig. [2] 
The data were recorded from March of 2011 to January 
of 2013. Also shown in Fig. [3] is the ground tempera- 
ture recorded during the same period. The data shows a 
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FIG. 1: (color online) Neutron daily counts percentage varia- 
tion (red) and atmospheric pressure (black) measured at the 
Sodankyla Geophysical Observatory in Oulu, Finland in 2012 
( |http: / /cosmicrays.oulu.n/| . The solid lines are two- week 
moving averages. 



significant drop of muon counts during the summer time 
while a higher counting rate seen in winter similar to the 
trend observed by the IceCube Neutrino Observatory 
and others (loj . 

It is therefore important to precisely understand how 
both the muon and neutron flux is attenuated in different 
layers of the atmosphere. Numerical simulations of muon 
and neutron flux variations at the surface of the earth 
have been carried out with varying air densities in the 
troposphere and stratosphere. The simulation software 
is developed based on the Geant4 package 
are reported in this letter. 



11|. Results 



In the simulation setup, a column of air 100 km in 
height and 50 km in diameter (composed of 70% Nitrogen 
and 30% Oxygen) is configured. The atmospheric air 
density is parameterized in the following equation (l2| . 
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FIG. 2: (color online) Muon daily counts percentage variation 
(red) and temperature (black) measured in Atlanta, Georgia 
from 2011 to 2013. Temperature data is from Atlanta Ful- 
ton weather station (http://www.wunderground.com/). The 
muon counts were recorded with a liquid scintillator detector 
installed on campus at Georgia State University. The solid 
lines are two-week moving averages. 
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FIG. 3: (color online) Simulated percentage variation of 
ground level neutron and muon counts as function of the per- 
centage variation of the troposphere density (black). Also 
shown in this figure is the daily neutron flux percentage vari- 
ation as a function of the atmospheric pressure from Oulu 
(red). 



where p is the air density in kg/m 3 , P is the pressure in 
kPa, and T is the temperature in Celsius. 

The primary cosmic particles impinging on the top 
of the earth's atmosphere consist of 79% protons and a 
small fraction of alpha particles and heavier nuclei 13 1 . 
In the present work, only primary protons are included 
in the simulation. The protons are launched vertically 
downward at the top of the air column with an energy 
distribution as described in [13j |. 



Based on the data provided by NASA [14|, the air den- 
sity in the stratosphere is larger in summer time and 
smaller in winter, which is opposite to the air density 
variation in the troposphere. This air density variation 
is modeled in our simulation by scaling the air density to 
match the seasonal variation in order to study its effect 
on neutron as well as muon flux changes measured at the 
surface of the earth. 

The empirical atmospheric data [13] shows a ±2.5% 
maximum variation of the troposphere density between 
winter and summer. This variation in density is modeled 
in our simulation by scaling the troposphere density from 
98% to 102% relative to the yearly average while keep- 
ing the same stratosphere density. Figure [3] shows the 
simulated neutron counts (black solid squares) variation 
at the ground level as a function of the percentage vari- 
ation of the troposphere density. Also shown in Fig. |3] is 
the measured neutron counts variation (red circles) from 
Oulu observation in 2012. The simulated result shows a 
very consistent trend with the data. This indicates that 
the neutron flux variation is primarily influenced by the 
modulation of the air density in the troposphere. One 
should not expect a perfect correlation between the mea- 
sured data and the simulated results for two reasons: (1) 



the simulation does not include the Oulu detector accep- 
tance and efficiency; (2) a constant primary cosmic ray 
flux is used in the simulation. 

While there is a strong anti-correlation between the 
troposphere air density and the neutron flux, there is lit- 
tle effect on the muon flux from the troposphere air den- 
sity variations as shown in Fig. [3] in black triangle points 
(except a noticeable effect on muon flux changes at the 
higher troposphere air density range). This is consistent 
with the results reported in [15j and led us to believe that 
the stratospheric air density may play a significant role 
for modulating the muon flux variation as seen in Fig. [5] 

In order to make a realistic simulation study, we 
modeled the stratosphere density variation according 
to the data from the Peachtree City Observation Sta- 
tion (Department of Atmospheric Science, University of 
Wyoming). In our simulation, stratosphere density is 
varied from 90% to 110% relative to its average while 
the troposphere density is kept constant. The results are 
shown in Fig. 2) A significant reduction of muon flux 
is clearly seen with increasing stratospheric air density, 
which is very consistent with our measurements at GSU. 
Based on the results in Fig. [3] and Fig. 01 it is clear that 
the effect of density fluctuation in the stratosphere re- 
gion is dominant on the muon flux variation comparing 
the modulation effect from the troposphere region. 

The modulation of muon flux from the variable strato- 
spheric air density can be further understood by study- 
ing the cosmic ray shower maxima altitude distributions. 
Primary cosmic ray showers mainly occur where the at- 
mospheric pressure is between 100 and 250 hPa, imply- 
ing that there is a certain density threshold (~0.2 g/m 3 ) 
necessary for maximum likelihood of primary cosmic ray 
interactions. Figure [S] shows the altitudes of the simu- 
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FIG. 4: (color online) Ground level muon flux percentage 
variation as a function of percentage variation of the strato- 
spheric air density. The red-circle data points are from GSU 
measurement and the black-square from the simulated results. 
We only sample muon particles with kinetic energy > 1 GeV 
at ground level which corresponds to the detector threshold. 
Note that the error bars in the figure are statistical only. 
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FIG. 5: (color online) Simulated cosmic ray shower maximum 
altitude (black square) as a function of the percentage varia- 
tion of the stratosphere air density. The red circle data points 
represent the mean altitude of the pressure range of 100 to 
250 hPa, which are taken from the Peachtree City Observa- 
tion Station (2011 - 2013). 



lated shower maxima (black squares) as a function of the 
stratospheric densities. Also shown in Fig. [5] is the vari- 
ation of mean altitude (red circles) of the stratospheric 
pressure range of 100 to 250 hPa extrapolated from the 
Peachtree City Observation Station data from 2011 to 
2013. Note that the larger error bars at low densities 
are due to low statistics. A remarkably good agreement 
is seen between the simulated shower maximum altitude 
distribution and the extrapolated pressure threshold re- 
gion. 

It is important to point out that the stratospheric air 



density variation indirectly reflects the seasonal temper- 
ature variation in the stratospheric region. The results 
shown in Fig. 0] suggest that muon flux is primarily influ- 
enced by the seasonal temperature variation in the strato- 
sphere. 

It is known that the temperature in the troposphere 
can fluctuate considerably within a day while the strato- 
spheric temperature only varies seasonally unless there is 
a sudden stratospheric warming 16]. The primary cos- 
mic ray particles mainly interact with the stratospheric 
nuclei and generate secondary cosmic ray particles at an 
altitude between 12 and 15 km. The mesons produced in 
these cosmic showers can either interact with the atmo- 
sphere or decay into muons. During the summer time due 
to the expanding atmosphere cosmic ray showers occur 
at higher altitudes. This means that the muons travel 
further to reach the surface of the earth and are more 
likely to decay leading to a lower muon rate in summer 
and a higher rate during winter. 

The simulated results in Fig. |4] are obtained by only 
varying the stratospheric density and keeping a constant 
tropospheric density. However, empirical atmospheric 
data [l4j show that the stratospheric and tropospheric 
densities vary inversely during summer and winter. The 
measured muon rate is higher compared to the simulated 
muon rate in summer times as a result of less modulation 
by the low troposphere density. This trend is also consis- 
tent with our muon simulation results shown in Fig. m 

In summary, we report in this letter for the first time 
the simulation studies of the effects of the tropospheric 
and stratospheric density variations on the cosmic ray 
muon and neutron flux reaching to the surface of the 
earth. Our simulation results show that the density vari- 
ations in the troposphere mainly influence the neutron 
flux while its influence on the muon flux is relatively in- 
significant. On the other hand, the stratospheric density 
dominates the muon flux which is in remarkably good 
agreement with observed seasonal variation of muon flux. 
It is, therefore, very important to have a longterm mon- 
itoring of both muon and neutron flux simultaneously at 
a global scale in order to study the dynamical change of 
the atmosphere. 
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